Matter density perturbations in modified gravity models with arbitrary coupling 

between matter and geometry. 
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We consider theories with an arbitrary coupling between matter and gravity and obtain the 
perturbation equation of matter on subhorizon scales. Also, we derive the effective gravitational 
constant Gef f and two parameters E and 77, which along with the perturbation equation of the matter 
density are useful to constrain the theory from growth factor and weak lensing observations. Finally, 
we use a completely solvable toy model which exhibits nontrivial phenomenology to investigate 
specific features of the theory. We obtain the analytic solution of the modified Friedmann equation 
for the scale factor a in terms of time t and use the age of the oldest star clusters and the primordial 
nucleosynthesis bounds in order to constrain the parameters of our toy model. 

PACS numbers: 04.50.Kd,95.30.Sf,98.80.-k 



I. INTRODUCTION 

There is accumulating observational evidence based 
mainly on type la supernovae standard candles [l| that 
the universe has entered a phase of accelerating expan- 
sion at a recent cosmological time scale. This expansion 
implies the existence of a repulsive factor on cosmological 
scales which counterbalances the attractive gravitational 
properties of matter on these scales. There have been 
several theoretical approaches Q , 0] towards the under- 
standing of the origin of this factor. The simplest such 
approach assumes the existence of a positive cosmological 
constant which is small enough to have started dominat- 
ing the universe at recent times. This model provides 
an excellent fit to the cosmological observational data 
Q and has the additional bonus of simplicity and a sin- 
gle free parameter. Despite its simplicity and good fit 
to the data, this model fails to explain why the cosmo- 
logical constant is so unnaturally small as to come to 
dominate the universe at recent cosmological times, a 
problem known as the coincidence problem and there are 
specific cosmological observations which differ from its 
predictions 

In an effort to address this problem two classes of mod- 
els have been proposed: The first class assumes that 
general relativity (GR) is a valid theory on cosmologi- 
cal scales and attributes the accelerating expansion to a 
dark energy component which has repulsive gravitational 
properties due to its negative pressure. The role of dark 
energy is usually played by a minimally coupled to grav- 
ity scalar field called quintessence\J\. Alternatively, the 
role of dark energy can be played by various perfect flu- 
ids (e.g., Chaplygin gas [Sj]), topological defects holo- 
graphic dark energy [1^, etc. The second class of models 
attributes the accelerating expansion to a modification of 
general relativity on cosmological scales which converts 
gravity to a repulsive interaction at late times and on cos- 
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mological scales. Examples of this class of models include 
scalar-tensor theories [llllT^. p^, f (R) modified gravity 
theories[14], braneworld models [1^, etc. 

Of these examples, f{R) theories have received much 
attention mainly due to the fact that they can provide a 
more natural explanation of the accelerating expansion 
of the Universe. At the same time many of these are able 
to satisfy the solar system constraints, see, for example, 
Ref. , and this is done without the introduction of ex- 
tra scalar field degrees of freedom, for a review of viable 
f{R) gravity models see [i2i],tl8i], while for the advan- 
tages and the drawbacks of this class of theories see [l9l | 
and [l^]. Also, the f{R) theories arise in a wide range of 
different frameworks: In quantum field theories in curved 
spacetime(2l|. in the low energy limit of the D = 10 su- 
perstring theory [Hi, in the vacuum action for the grand 
unified theories, etc. 

Recently another approach was proposed, first in 
[13 1 01 ^"^^ later in [2^ , namely, to add as a nonminimal 
coupling to matter a general function of the Ricci scalar. 
Such nonlinear couplings of matter were analyzed in the 
past |26i in the context of the accelerated expansion of 
the Universe and in Refs. [l^l, [111 it was shown that 
these theories do not correspond to scalar-tensor grav- 
ity and therefore have a yet unexplored and non-trivial 
phenomenology. Furthermore, such theories due to the 
explicit nonminimal coupling between the matter content 
and gravity have a nonvanishing covariant divergence of 
the matter energy-momentum tensor T^^ [2^, [SOI, [3l| . 
which means that there is an energy exchange between 
matter and gravity. 

In Sec. II we will derive the background equations 
of motion for a general Lagrangian density \fi{R) + 
G{Lm)f2{R) and we will discuss some implications about 
the possible choices for the matter Lagrangian Lm- Hav- 
ing more than one Lagrangians to describe the same fluid, 
e.g. Lm — p or Lm — — p for pressureless dust, is not a 
problem in GR as all choices will give the same equa- 
tions of motion, but in theories with a nonminimal cou- 
pling, different choices in general will give different phe- 
nomenology. Nevertheless, we will choose a Lagrangian 
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that depends only on the density Lm = Lm{pm) in or- 
der to keep the analysis as simple as possible, but quite 
general at the same time. 

In Sec. Ill we will derive the matter density perturba- 
tion equations for the case that the matter Lagrangian is 
a general function of the density p. Being able to predict 
the evolution of the matter perturbations for such mod- 
ified gravity models is necessary in order to compare the 
theory against the linear growth data. Such data are in 
principle able to differentiate between GR and modified 
gravity, with the main reasons being the evolution with 
time of the gravitational constant 32 1, [ssj . 

We will also evaluate two parameters rj and S — o (l -I- 
77/2) for this particular model, following Refs js^.js^ 
where the first quantifies the strength of an anisotropic 
stress while in the latter, q quantifies the deviation of the 
gravitational constant today. These parameters can be 
used to constrain the theory with weak lensing observa- 
tions and if such surveys measure deviations from GR, 
for which iJl,r])\GR = (1,0) then those two parameters 
could be used to differentiate between modified gravity 
and dark energy models. 

Finally, in Sec. IV we will use a toy model with La- 
grangian density R/2 -I- (l-f XR)G{L,n), where G{Lm) is 
given by G{Lm) — AL^^ and L„i = —p- This model 
has the energy-momentum tensor of pressureless dust, 
but exhibits nontrivial phenomenology compared to its 
counterpart from GR. Also, we will present an analytical 
solution to the modified Friedmann equation and give the 
scale factor a in terms of the cosmic time t. In order to 
provide some bounds on the parameters A and n of our 
toy model, we will use the following two observations: the 
age of the oldest star clusters and the Big Bang Nucle- 
osynthesis (BBN) bounds on the gravitational constant 
Gn, which arc of the order of 10%. 



II. BACKGROUND EQUATIONS 



The action we will consider is 



S 



d X\J—g 



h{R) + G{L„,)f2iR) 



(2.1) 



where g is the determinant of the metric g^^i, and fi — 
fi {R) , i = 1,2 are arbitrary functions of the Ricci scalar 
i?, while G{Lm) is an arbitrary function of the matter 
Lagrangian density L^. We use the metric signature 
(-,+,+,+). 

Varying the action with respect to the metric g^j/ we 
obtain the field equations as 

[Fi(i?) + 2G{L„,)F2{R)] R^, - ^fi{R)g^. + 
+ [5^,n - V^V,] (FiiR) + 2G{L„,)F2{R)) - 



f2{R)K{L,^)T^ 



(2.2) 



where Fi = dfi/dR, i = 1,2, i?^,^ is the Ricci tensor and 
KiL^) = dG{Lm)/dL^. When /i(i?) = /(i?), ^(i?) = 



1 and G{Lm) = Lm then we recover the field equations 
for /(i?) gravity. 

Finally, we have defined the energy-momentum tensor 
of the matter as 



T, 



2 Sj ^Lrn) 
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(2.3) 



Assuming that the matter Lagrangian density Lm de- 
pends only on the metric and not on it's components, 
we obtain T^^ — Lmgi^v — 2dLm/dg'^'^. In what fol- 
lows we will assume that the matter Lagrangian is an 
arbitrary function of the energy density p only, so that 
Lm = Lm (p) ■ Then the energy-momentum tensor can be 
written as 



T^" = -p^u'^u'' + I L 
dp 



dp 



(2.4) 



where the four- velocity — dx^ /ds satisfies the con- 



dition 



= — 1. In deriving Ea. (|2.4p we have used 



Eq. ()2.3|) . the relation Sp = \p{gfiv + u^Uu)5g^'^ and 
the conservation of the matter current W^ipu") — 0. 
Note that the difference in the sign of the first term be- 
tween our Eq. (|2.4p and the corresponding one of Ref . [3l| 
is due to the difference in the normalization condition, 
u^Ufj_ — —1 for us. 

In general the Lagrangian of a perfect fluid (on shell) 
is equal to the pressure L,„ — p, ie for dust it is equal to 
zero. However, this is not the only choice as the La- 
grangian Lm = —p also reproduces the same energy- 
momentum tensor T^^ = pu^u^. This is not a prob- 
lem in GR as both choices give the same equations of 
motion, but in theories with a nonminimal coupling, dif- 
ferent choices give different phenomenology. This might 
be a problem since our choice for the Lagrangian be- 
comes physically meaningful without having any criteria 
for choosing it, besides the observations of course. 

This fact has created some debate on the literature 
lately [23|,[29|,[33| regarding which is is the best choice 
for the Lagrangian Lm for the perfect fluid. Nevertheless, 
we choose a Lagrangian that depends only on density 
Lm = Lm{pm) in Order to keep the analysis as simple as 
possible but quite general at the same time. 

In a flat Friedmann-Robertson- Walker metric with a 
scale factor a(i), we obtain the background (zero-order) 
equations: 



^F^+2GF2\H^ = f2Kpm + ^[{Fi+2GF2)R- h] 



-f2{G-KLm) - 3H-{Fi+2GF2) 
dt 



2 [i^i -I- 2GF2] H = f2Kpm 

2GF2] 



(2.5) 



(2.6) 



where H = R = 6{2H^ + H) and a dot represents a 
derivative with respect to the cosmic time t. 
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In order to get the evolution equation for the ideal fluid 
we take the covariant divergence of Eq. (|2.2p and by 
using the purely geometrical identities ^91: V^G^i, = 
and {nW^ - V^n)F = R^^WF we get|311 



For the case at hand (|2.7p gives 



(2.7) 



Prn 



-3K[f2Hprn + 2F2{L^ + p„r){AHH + H)] 

h{K + [L^ + Prn)K'{Lrn)L'^iPrn)) 



(2.8) 



where the prime denotes differentiation with respect to 
the particular argument, i.e., K'{L„i) = dK/dL„i while 
L'miPm) = dL„Jdp„i. Obviously, when G = and 
f2{R) — 1 we recover the usual energy conservation equa- 
tion of GR and f{R) theories. 



III. PERTURBATION EQUATIONS 

We will consider the following perturbed metric with 
scalar metric perturbations $ and ^' in a longitudinal 
gauge: 

ds^ = -(1 + 2<^)dt^ + a{tf{l - 2^)5rjdx'dx^ (3.1) 

The energy-momentum tensor of the nonrelativistic mat- 
ter is decomposed as Tq = —{pm + 6pm) and T° — 
—PmVm,a, whcrc Vm is & vclocity potential. The Fourier 
transformed perturbation equations are given in the Ap- 
pendix by Eqs. (jM]), and 

Next, we define the gauge invariant matter density per- 
turbation 5m as 



mv 



where 



avr, 



(3.2) 



(3.3) 



Following the approach of Ref. [35| , we use a subhori- 
zon approximation under which the leading terms cor- 
respond to those containing and 5pm- Terms that 
are of the form or $ are considered negligible rel- 

ative to terms like (fc^/a^)<i> for modes well inside the 
Hubble radius (fc^ ^ a?H^). Under this approximation 
Eqs. ((X3|) . (|M)) and ((^ yield a second order difl'erential 
equation for (5, given in the Appendix by Eq. (|A5|) . 

Equation (|A5[) compared to the corresponding one of 
GR or f{R) theories (see for example [11], [s^, [3^) has 
an extra term which is a function of the matter content 
of the theory. This fact, as we will see in what follows, 
will force Newton's constant Ge// to be dependent on 
the matter content, something which is in sharp contrast 
with f{R) theories, but not scalar-tensor gravity where 
Geff may depend on the scalar field (j) through it's non- 
minimal coupling F{(j)). 



Next, we write 5Fi and 5F2 as 

5F, = F,,R SR 



(3.4) 



where i = 1,2, the comma denotes derivative with re- 
spect to the Ricci scalar R and SR^ under the subhorizon 
approximation, is given by: 



(5i?~-2^(2^'-$) 



(3.5) 



Making these substitutions and using the subhorizon 
approximation in Eqs. (|Aip and (|A2p we get a set of 
equations for <i> and ^ given by Eqs. (|A6|) and (jA7[) . 
The next step is to express $ and in terms of 6m ■ This 
can be done by solving the system (|A6[) and (|A7p for $ 
and Doing so we find <i> and ^' under the subhorizon 
approximation, given by Eqs. (|A8[) and (jA9p . 

From the equation for $ (|A9|1 we can define a Pois- 
son equation in the Fourier space and attribute the extra 
terms that appear on the right-hand side to an effective 
gravitational constant Geff - Doing so, we get the gravi- 
tational potential 



-4:TTGeff-^SmPn 



(3.6) 



51) can 



where Geff is given by (jAlOp . 

Hence, the equation of matter perturbations 
be written compactly as 



Sm + 2HSm - MGeff + P)Pm5m - (3.7) 



where the parameter [3, after using Eqs. (|A8|) and (jA9[) . 
is given by (jAll[) . As it can be seen from Eq. (jAll[) in 
the minimal coupling case [f2iR) = 1 and G{Lm) = Lm] 
we get the expected result /? = 0. 

We also define two parameters 77, which characterizes 
the strength of the anisotropic stress 



(3.8) 



and q via = —inGoqpmSm, where Go is the value of 
the gravitational constant measured by the solar system 
experiments today. Then these two parameters are given 
by Eqs (^12)) and (^131) . 

In order to confront the modified gravity models with 
weak lensing observations it is necessary to calculate the 
potential ^wl that characterizes the deviation of light 
rays. Under the assumptions that photons follow null 
geodesies and that the geodesic deviation equation holds, 
the lensing potential is ^wl — $ -I- 5*; see Ref. [39] for 
a quick derivation and a discussion. However, in theo- 
ries with nonminimal coupling to matter this might not 
always be the case. In Ref. [27] it was shown that for the- 
ories with Lagrangian /i(i?)/2 + (1 + Xf2{R)) Lm, non- 
minimal coupling corrections to the Maxwell equations 
which would otherwise affect standard optics, can actu- 
ally affect only long, compared to the radius of the cur- 
vature of spacetime, wavelengths. When this is not the 
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case, i.e. we consider the high-frequency hmit, photons 
are transverse and propagate along null geodesies. 

Having this in mind, we define a combination of pa- 
rameters I] = (7(1 + 77/2), for which 



WL 



$ + ^ 



-87rGo^Pm(5mS (3.9) 



where A, A and n are constants. From Eq. (|2.4p we 
see that the energy-momentum tensor then is given by 
T^v — pUfj,Ui, as usual. The modified Friedmann and 
continuity equations can easily be found from Eqs. (j2.5p 
and (j2.8[) . For the continuity equation we find 



3n{a) + an'{a) = 



(4.1) 



where the explicit form of S is given by Eq. (|A14p . 
Again, in the case of the minimal coupling [f2{R) = 1 
and G{L„i) — L„i] we get the expected result for S: 



1 



SnFGn 



(3.10) 



which agrees with the result of Refs. [3^,[3J] for simple 
f{R) theories. We should remind the reader here that 
the formulas derived for the weak lensing potential have 
been derived under the assumption that light travels in 
null geodesies, which may be the case only for certain 
models and in the high-frequency limit. 

Therefore, this class of theories may have non-trivial 
effects on the weak lensing that could potentially be ob- 
servable. However, this would also require one to work 
out the most general case for the relation between the 
lensing potential and the metric perturbations. Since 
this is beyond the scope of this analysis, we have left the 
most general case for future work. 

In order to recover f{R) gravity for all the results 
we have mentioned, we only have to set fi{R) = fiR), 
f2{R) = 1, and G{Lm) = This means that F2{R) — 
0, K{L,n) = 1, and K'{L„i) — 0. Then, for example, Eq. 
(|XTOl) yields 



1 1 



4-|i5m 



SttF 1 



3-^m 



(3.11) 



where 



RFi 



m 



also the anisotropic parameter 77 defined in (j3.8p and 
given by (|A12p becomes 



+ 2^m 



(3.12) 



with both equations, (jS.lip and (|3.12l) . being in agree- 
ment with the standard results from f{R) gravity [35|. 



IV. A TOY MODEL 

As an example we will now consider the case where 
we have the Lagrangian of GR plus a specific nonmin- 
imal coupling of gravity to pressureless dust. Specifi- 
cally, for our toy model we will assume that fi{R) = R, 
f2{R) = 1 + \R, Lra{p) = -p and G(L„) = ^Li+", 



where the derivatives are written in terms of the scale 
factor a and we have defined p{a) = 3HQn{a). 

Equation (|4.ip can as usual be solved for fl{a) and with 
initial conditions ~ 1 we get 



n{a) 



1 



(4.2) 



For the modified Friedmann we find, after using 



(-3)"i/2«+2^ 



2(-l)"3"+M(3n + 2)XHl"+^ + a3"+3 



(4.3) 



By demanding that H{1) = Hq we can solve for A and 
after substituting back to the previous equation we find 



H{ar 



Hi 



a3"+3 (1 - QHli'in + 2)A) + 6H^{3n + 2)A 



(4.4) 



Obviously, when A = and n — we get the usual, 
from GR, Friedmann equation for a matter dominated 
universe. If we expand (|4.4p for A <C 1 we get 



6a-s("+i) (a3»+3 _ 1) H^l^n + 2)X + (A2)](4.5) 

Thus, when A ^ 1 such models can be thought of as a 
sum of ideal fluids with the first term being attributed 
to a matter fluid component with constant equation of 
state w — n, while the second to a dark energy fluid with 
variable equation of state w = w{z). 

We can use Eq. (|4.4p to solve for a{t) in terms of the 
cosmic time t with the initial condition a (to) = 1, where 
to is the present time (the age of the universe). Doing so 
we find a{t) in implicit form: 

^ ^ a(t)3"+3m3F2 (i 1, 1; 2, 2; -a{tf^+^m) 
\A^r+T(6n -h 6) 



l-6Hg(3n+2)A 



and 3F2 is a hypergeometric 



where m 

function. Then the age of the universe is 

mgi^a (i, 1,1; 2, 2; -to) 



t{a = l) = to 



6Ho^/m+l{n + 1) 



(4.7) 



Demanding that the universe should be at least old 
enough to accommodate the oldest star clusters [1^ , then 
it is necessary to have to > 12Gyr. Equation l|4.7p also 
implies that A and n should satisfy certain conditions in 
order to have a real value for Iq. Specifically, we find 
that these conditions are A > and ri>— |orA<0 
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combine the two constraints there is no allowed param- 
eter space at all for our model, there are actually val- 
ues for {n, A) that satisfy both observational constraints, 
e.g. (n. A) — (—2/3,0.015873), but they require extreme 
fine-tuning to achieve agreement with the observations 
and some of these values exhibit unphysical behavior for 
Geff at early times, such as singularities. Therefore we 
may conclude that our simple toy model is in some stress 
with observations as the allowed parameter space is very 
small and even then for some of the allowed parameter 
values, the model may exhibit unphysical behavior. This 
of course does not mean that the whole class of theories 
with a nonminimal coupling is problematic, as one may 
well choose more suitable forms of the functions fi{R), 
f2{R) and G{Lm) that satisfy both observational con- 
straints. 



CONCLUSIONS 



FIG. 1; The allowed parameter space of (n, A) such that to > 
l2Gyr, where A is in units of I/Hq. The embedded plot shows 
the region around A ^ 1 near the upper bound for n. 

and n < — I . In Fig. [T] we show the allowed parameter 
space of (n, A) such that to > 12Gyr. As it can be seen 
from the embedded plot there is an upper limit to n that 
corresponds to A = and to — , with the latter 

giving n < —0.256 after we use the bound to > l2Gyr. 

Next we will use the primordial nucleosynthesis (BBN) 
constraints on the variation of the gravitational constant 
to further constrain the parameters A and n. The effect of 
the variation of Gef / can be constrained from BBN to be 
of the order of 10%, see for example Ref. JdJ which gives 
= 1.09±[5;f|. The effective gravitational constant 
Geff can be calculated from (jAlOP by using (|4.2p and 
(|4.4p . If we set Geff{a = 1) = Gq then the gravitational 
constant is given by Eq. (jAlSp . 

First, we will consider the following two cases, for A = 
and n — 0. In the first case the gravitational constant is 
simply given by Geff = Goa~^" and by using the BBN 
constraint of Ref. [4l| we find -0.0018 < n < 0.0045 
while in the second case by using a series expansion for 
small A, which is to be expected from the Solar System 
tests, we find 

Geff /Go = l+(^2fcVij2 + ^'^^f/^° - 9) ij2A+0 (A^ 

(4.8) 

In this case using the BBN constraint we get —8.9 • 
10-29 < XH^ < 2.8 • 10-2«. In the general case, the 
allowed parameter space, shown in Fig. 2, is rather 
complicated. When n > then A is constrained to be 
-2.8 • 10^2^ < XH^ < -1.8 • 10^27^ ^YiWe when n< we 
cannot give a bound for A in closed form. 

It should be noted that while it seems that when we 



We derived the matter density perturbation equation 
and the effective gravitational "constant" Geff for the 
action (|2.ip . Our analysis covers f{R) models with an 
arbitrary coupling of the matter content to gravity for 
any matter Lagrangian that depends on the fluid density 
L,n = Lm{pm)- In Order to recover the ideal pressure- 
less fluid we chose L,n — —p- Having chosen some other 
dependence for the Lagrangian our results would be dif- 
ferent since the equations of motion depend explicitly on 
the choice of the Lagrangian. As we have mentioned ear- 
lier, having more than one Lagrangian to describe the 
same fluid, eg Lm = p or Lm = ~p for pressureless dust, 
is not a problem in GR as all choices will give the same 
equations of motion, but in theories with a nonminimal 
coupling, different choices in general will give different 
phenomenology. 

We also evaluated two parameters 77, which quantifies 
the strength of an anisotropic stress, and S = (7(1-1-77/2), 
where q quantifies the deviation of the gravitational con- 
stant today. These parameters can be used to con- 
strain the theory with weak lensing observations and 
if such surveys measure deviations from GR, for which 
(E,77)|Gfl; = (1;0) then those two parameters could be 
used to differentiate between modified gravity and dark 
energy models. 

Finally, we used a toy model with Lagrangian den- 
sity R/2 + (1 -h XR)G{Lm), where G{Lm) is given by 
G{Lm) = AL"^^ and = —p. This model has the 
energy-momentum tensor of pressureless dust, but ex- 
hibits nontrivial phenomenology compared to its coun- 
terpart from GR. We presented an analytical solution to 
the modified Fricdmann equation and gave the scale fac- 
tor a in terms of the cosmic time t. However, our toy 
model is in some stress with observations as the allowed 
parameter space is very small. In order to have an old 
enough universe, so as to accommodate the oldest ob- 
served clusters, the required values for A are much larger 
than those allowed by the BBN constraints (see Figs. 1 




and 2). However, this fact should be interpreted with 



care since Ref. [4l| used a simple modification of the 



Friedmann equation to derive the BBN constraints we 
used. A full analysis of the implications of modified grav- 
ity models with arbitrary coupling between matter and 
geometry during nucleosynthesis was beyond the scope 
of this analysis and is left for future work. 
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APPENDIX A: THE PERTURBATION 
EQUATIONS 
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In this Appendix we have gathered some of the equa- 
tions that if left in the main text they would make the 
analysis difhcult to follow. 

The Fourier transformed perturbation equations are 
given by 



-f2KSpm - HSFi + 2GSF2)H + 6F2KH{2Lm^ + 2p™$ - SpmLUPm)) 

-/2<Sp,„(i™ + p„,)K'{L^)L'JPr,r) - 3H\SFi + 2{Fi'l> + G{6F2 + 2^2$) + F2i^(-4(L„ + p,„)$ 
+Jp,„L;,(p„0))) - i^H^iSFi + 2(-Fi* + G(SF2 - 2F2*) + F2K{p„, - i,n($ - 2*) 

+dp„MraiPrnm - 3(A + 2GF2 + 2F2K L'„Xp,n) p'n.)-^ + 6F2X(L„ + p„,)* 

+3i/(-2$(Fi + 2GF2 + 2F2KL'„XPn.)p;n) + 2F25p,nL[^{pmfpmK'{L„,) + SF, + 2G6F2 - 
2(Fi + 2F2G)* + 2K{6F2LUPm) + 5 p^^iALLiPm) + F2p;„L:'„(p„)) 

+F2{L'^ipm)5p,n + {Lm + + 4*)))) = (Al) 



(j) _ \I/ 



SFi + 2G (Lm) - ^)F2 + 5F2) + 2K (L„0 F25p^L'^ jpm) 

Fi 



(A2) 
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-6p^ - SHSpm + p (^3* - — u„j = {~6a^H + AHH)F26f2(,L^ + Pm)K^ - f2{2F2{l2H^^Lm + Pm)a^ 

.... . /)3vl/ 

+3((i:,„ + Pm){2<^H + i7(3$ + 4^-) + — ) - H8p„,{L'„,{p^) + l))a'' - k'i^ - 2*)(i:,„ + p™) 
+i?(-12ij(5p,„(L;^(p„) + \)a^ + (3($ + 4*)a2 - 4fc2* + 2$(12ija2 + e))L,n + (3($ + 4*)a2 
-AkH + 2$(12ija2 + A:2))p„,)) - aHhiL^ + Prn))if' + fi{K{5p„,{L„, + p„,)K' {L^)L'^{p^) 

)L':n{Pm)))) 

-PrJpmiLm + p^)K' {Lmf L'^{pmf)) I {a^Kf^{K + 5 p„,K' {L„,)L'^{p„,))) (A3) 
+pf2PmVmK\L^)L'^{p„,)a^ + {K{L^Sf2 + (!>a{H^ + 3iji7 + H)F2p^v^) 

+ f2LmSpmK'{Lra)L'^iPm))a^)/ {o" f2{K + Sp^K' {Lm)L'^iPm))) (A4) 

I 



Under the subhorizon approximation the differential 
I 



Using the subhorizon approximation in Eqs. (|Af p and 
(|A2p , we get a coupled set of equations for $ and 5* : 

\{~ l25mPni{K + (L„i + Pm)K' {Lm)L'^iPm))a'^ 

-2fc2(*Fi + F2(2G* + - 2W)Lm 

+Pm{^~2^-5„X„APrn)))))a^ 

+2(Fi^fl + 2F2^flG)fc4($-2*)) = (A6) 



((($ _ + 2F2(G($ - + /f(5„,/9,„£:„(p,„)))a' 
+2(i^i^7j + 2i^2.flG)fc2($ - 2^))/{a^Fi) = (AT) 

By solving the system of equations (jA6p and (|A7p for $ 
and 5* we find 

2 

^ = -^S,nP,nif2iFl + 2GF2){K + (i„ + p^) 

K'{L,r,)L'^{prn))a^ + 2k^-KF2{F, + 2F2(G 
+X(i™ + p,„)))L:„(p„) + (^^1,^ + 2F2,flG)/2 

(if + (L„, + p„0^'(im)i:n(Pm))))/ 

(2(i^i + 2GF2)((Fi + F2(2G - K{L^ + p,n)))a^ 
+3iFi,R + 2F2.,RG)k^)) (A8) 



equation for the growth of perturbations S is 



(A5) 



and 

* = -p'5„,p,„(/2(i^i + 2GF2){K + (L,„ + p,„) 

if'(L„0L:„(p„0)a2 + 2e{KF2{Fi + 2F2(G - 
2if(L„ + p™)))i;,(p„) + 2F^MK + [L 
K'{L^)L'^{p^)) + AF2,RGf2{K + (L„ + p„) 
K'{L,r,)L'^{p„r))))/{2{Fi+2GF2){{Fi+F2{2G^ 
K{Lm + Pm)))a^ + HF^M + 2F2MG)k^)) (A9) 

By using (jA9[) to define a Poisson equation in the 
Fourier space we get Ge/ / 

Geff = ^(/2(-Fl + 2GF2)(if + (Lrn + Pm) 

K'{Lm)L'^{pm))a^ + 2k\KF2{Fi + 2F2(G - 
2K{Lm + Pm)))L'^{Pm) + 2FiMK + (i™ + Pm) 

K'{Lm)LlJpm)) + AF2,RGf2{K + (i„ + Pm) 
K'{Lm)KniPm))))/{{Fi + 2GF2){{F^ + F2(2G - 

K{Lm + Pm)))a' + 'S{Fi^R + 2F2MG)k')) (AlO) 
The parameter [3, which corresponds to the last term 



5m + 2H5m H 5 ^{{Lm{K - 5mPmK'{Lm)L'„i{fim)){f25rnPmK'{Lm)L'^{Pm)a^ 

+2eK{<^> - 2^)F2))/{a^K^f2Pm)) ^ 

I 
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of Eq. 



/3 



5j), after using Eqs 

1 e 



8j) and is given by 



A-K a' 



Lrn{f2{F2iK^ + 2GK' {L^)L[^{p^))a' 



+ {Fia^ + 3{FiM + 2F2,BG)k^)K'{L^)L'„^{p,n)) 
-&eK^FiL'^{pm)))/{Khp„^{~F,a' 
+F2{-2G + KL,n + Kp^)a^ 

-3(Fi,fl + 2F2MG)k^)) (All) 



The two parameters 77 and q defined in the main text 
are given by 

77 = 2^(2XF2(Fi + F2(2G - KL„, - Kp„,))L'„M 
+ iFi,B. + 2F2,RG)h{K + L,^K'{L,r,)L'^{p^) 

+p^K'{L,r,)L'^{pra)))/{f2{Fl + 2GF2){K + (i,„ + pm) 

K'iL„,)LUpra)) + 2fcVa'((Fi,i^ + 2F2,RG)f2{K 

+LmK' {Lm)L[^{pm) + PmK' {L„i)L'^{pm)) 

-KF2(Fi + 2F2{G + XL™ + K p„,))L'^{p„,))) 



(A12) 



and 



1 



9 = 



-(/2(Fi+2GL^2)(i^ + (L„ 



SttGo 

K'{L,n)L'^{pm))a^ + 2e{-KF2{Fi + 2F2(G 

+ P„)))L;„(p„) + Fi,Rf2{K + {Lrn + Pm) 

K'{L,^)L[^{p,^)) + 2F2MGh{K + (L„, + p„) 
if'(L„)L:„(p„))))/((i^i + 2GF2)((Fi + F2(2G - 
K{L„, + pra)))a'' + i{Fi^B, + 2f^2,flG)fc2)) (A13) 



The expHcit form of E, used to confront the modified 
gravity models with weak lensing observations, is given 

by 



S = {h{Fl + 2GF2){K + [L^ + Prn)K'{Lrn)L'^{Prn))a^ 
+ 3k^iiFi,R + 2F2,RG)Kf2 - {Lra + p^m) 

{2K^Fl - [F^.R + 2F2.RG)hK'{F^))L'^{pm)))/ 
{8Go7t{Fi + 2GF2){{Fi + F2(2G - K{L„, + p^)))a^ 
+3(i^i,fl + 2F2,RG)e)) (A14) 



In the case of our toy model of Sec. IV, Ge// can be 
calculated by using Eq. (jA10[) and is given by 



Geffia) = Go(a(18(n + 1)XH^ - l)(-3A(9a2(n + 1) - 8k^{3n + 2)\)H^{6(3n + 2)XH^ - l)a^("+i) + 
(a^ - 2A:2a)(1 ~ 6(3n + 2)\H^^)^a^''+^ + 36(3n + 2)X^{a^{3n + 4) - 2k^{3n + 2)X)H^))/ 
(((6(3n + 2)XHl - l)a3"+3 + 6Aij2)(a3"+3(6(3n + 2)XHI - 1) ~ 6(3n + 2)XHlf 

(1 + A(-— 2 + 54(n + l)(3n + 2)XHI - 9n + i)Hl)) (A15) 
^0 
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